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Abstract. Using the test-field method for nearly irrotational turbulence driven by spherical 
expansion waves it is shown that the turbulent magnetic diffusivity increases with magnetic 
Reynolds numbers. Its value levels off at several times the rms velocity of the turbulence multi- 

' plied by the typical radius of the expansion waves. This result is discussed in the context of the 

, galactic mean-field dynamo. 
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The galactic dynamo is believed to be powered by supernova-driven turbulence. This 
type of forcing does not directly produce vorticity; it can only be produced indirectly 
through oblique shocks, i.e. through the baroclinic term. The aim of this work is to assess 
whether vorticity is actually important for the dynamo. 

The galactic magnetic field has a strong large-scale component which is generally 
believed to be due to a mean-field dynamo of aQ type that is governed by the equation 

^ 

— = Vx(UxB + £- ?y^o J), where £ t = n,,77, - //,,//,,7, (0.1) 

is the mean electromotive force, J = V x B/fM> is the mean current density, U is the 
mean flow, and fiQ is the vacuum permeability. In order to assess its efficiency, one needs 
to determine the tensors ctij (the "a effect") and rjij (turbulent magnetic diffusivity). 
Note that aiij is a pseudo tensor and non-vanishing diagonal components can only be 
constructed from a combination of polar and axial vectors, and would therefore be van- 
ishing in the absence of stratification and rotation. The rjij tensor, on the other hand, 
■ does not require this, and it should be finite even in the completely homogeneous case. 
\ This is the case considered in the present study, which is a necessary intermediate step. 
For homogeneous flows rjij is an isotropic tensor, which we write as rjij = rj t 5ij, where 
rjt is the turbulent magnetic diffusivity. A relevant concern in mean-field theory is that 
turbulent transport coefficients such as rjt must stay finite even in the limit of large 
values of the magnetic Reynolds number, defined here as R m — Urms/rjkt, where w rm s is 
the rms velocity of the turbulence, and fef is the wavenumber corresponding to the scale 
of the energy-carrying motions. Given the importance of a possible R m dependence, it 
is necessary to perform so-called direct simulations, where no subgrid scale modeling is 
used. This implies that we must make compromises regarding the strength of the forcing 



> 



ON 

o 



X 



and consider only subsonic flows. Following earlier work of Mee & Brandenburg (2006) 
we consider a flow driven by random expansion waves of radius R = 2/kf (not to be 
confused with the magnetic Reynolds number R m ) and determine rj t using the test-field 



method of |Schrinner et al. (2005) in the implementation of Brandenburg (2005) 



The evolution of internal energy and hence entropy is not relevant to our question 
about turbulent transport coefficients. Therefore we consider an isothermal equation of 
state where the pressure p is proportional to the density p with p = pc^ , with c s being 
the isothermal speed of sound. We adopt a Gaussian potential forcing function / of the 
form f(x,t) — V</> with 4>(x,t) — iVexp{[a; — Xf (i)] 2 /i? 2 }, where x = (x,y,z) is the 
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Figure 1. Dependence of rjt on P m = 1. 

position vector, ccf(t) is the random forcing position, R is the radius of the Gaussian, and 
N is a normalization factor. We consider a time dependence of Xf with a forcing time 
8t[ OTCC ~ (urmgfcf)" 1 that defines the interval during which Xf remains constant. We use 
the Pencil Code (http://pencil-code.googlecode.com) which is a non-conservative, 
high-order, finite-difference code (sixth order in space and third order in time) for solving 
the compressible hydrodynamic and hydromagnetic equations. 



In Figure[T]we plot the dependence of rj t on R m . Following earlier work of Sur et al. (2008) 
we normalize r/ t by r/ t o = it rms /3fcf. Note that, for low values of R m , rj t increases propor- 
tional to -R™ with n between 1/2 and 1. For larger value of R m , rjt seems to levels off 
at a value of about 20 times rjto ■ Expressing this in terms of it rms and the typical radius 
R of the expansion waves, we find that rjt 4it rms i?. Note also that rjt is always posi- 
tive, in contrast to analytic predictions for irrotational turbulence using the first-order 



smoothing approximation (Krause & Radler 1980) 



Based on these results we can conclude that nearly irrotational turbulence is at least as 
efficient as vortical turbulence in diffusing mean magnetic field. Clearly, our study is still 
at a preliminary stage. It is important to clarify a possible dependence of our results on 
the microscopic magnetic Prandtl number and, in the nonlinear regime, on the magnetic 
field strength. Next, we need to consider the case with rotation and stratification which 
should then lead to an a effect, as well as turbulent pumping. This would provide an 



opportunity to compare with early predictions by Ferriere (1992) for this type of flows 



Acknowledgements 

We acknowledge the allocation of computing resources provided by the Swedish Na- 
tional Allocations Committee at the Center for Parallel Computers at the Royal Institute 
of Technology in Stockholm and the National Supercomputer Centers in Linkbping. This 
work was supported in part by the European Research Council under the AstroDyn 
Research Project 227952 and the Swedish Research Council grant 621-2007-4064. 



References 

Brandenburg, A. 2005, Astron. Nachr., 326, 787 
Ferriere, K. 1992, ApJ, 389, 286 

Krause, F., & Radler, K.-H. 1980, Mean-field magnetohydrodynamics and dynamo theory (Perg- 

amon Press, Oxford) 
Mee, A. J., & Brandenburg, A. 2006, MNRAS, 370, 415 

Schrinner, M., Radler, K.-H., Schmitt, D., et al. 2005, Astron. Nachr., 326, 245 
Sur, S., Brandenburg, A., & Subramanian, K. 2008, MNRAS, 385, L15 



